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ABSTRACT 

A  system  is  describe!  which  incorporated  a  Hewlett-Packard 
85  "desk-toe"  computer  to  control  a  frequency  synthesizer 
and  read  the  output  of  a  l^ck-in  analyzer  to  measure, 
display  and  record  tha  resonant  frequencies,  amplitudes,  and 
quality  factors  of  several  modes  of  an  acoustical  resonator. 
The  system  is  capable  of  locating,  measuring,  and  tracking 
the  resonant  medes  as  parameters  which  affect  sound  speed 
and  attenuation  are  varied.  An  algorithm  for  rapidly 
fitting  "good  quality"  measured  data  to  a  resonance 
lineshape  is  described  which  determines  quality  factors  to 
precisions  cf  better  than  9.1  percent,  amplitudes  to  better 
than  0.01  percent  ind  center  frequencies  to  better  than 
O.lppm.  Sample  output  is  prDvided  for  the  lowest  three 
plane  wave  modes  of  a n  air  filled  cylindrical  resonator  in 
the  temperature  ranga  of  -15  to  25  degrees  Celsius. 
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A.   BACKGROUND 

Many  applications  in  physics  and  engineering  require  the 
ability  to  determine  tha  cantar  frequency  (fo) ,  maximum 
amplitude  (A),  and  quality  factor  (Q)  of  a  resonance,  A  few 
examples  of  this  need  arise  ia  the  measurement  of  the  speed 
of  sound,  in  the  reciprocity  calibration  of  electroacoustio 
transducers  [  Bef .  1],  and  in  the  measurement  of  sound 
absorption. 

Prior  to  the  ai/ ent  of  computers  and  instrumentation 
interfacing,  sever?  liaitations  were  imposed  upon  the 
experimental  precision  obtainable.  The  variation  of 
amplitude  at  resonanre  is  seconl  order  in  frequency  so  that 
a  small  uncertainty  in  the  maximum  amplitude  results  in  a 
greater  uncertainty  in  the  frequency  corresponding  to  the 
maximum  amplitude.  Uso,  the  3  is  calculated  from  the  ratio 
of  a  large  number  (renter  freqiency)  to  the  difference  of 
two  large  numbers  whirh  may  differ  by  a  snail  amount  (the  3 
dB  down  frequencies),  resulting  in  an  uncertainty  in  Q  which 
can  be  quite  large. 
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This  project  was  designed  to  produce  an  automated,  high 
precision  system  capable  of  obtaining  and  tracking  the 
center  frequency,  maximum  amplitude  and  2  of  several 
acoustical  resonant  modes  in  any  given  resonator  as  some 
external  parameter  affecting  these  values  is  varied.  The 
approach  was  broken  up  into  four  task  groups. 

1.   Acoustical  Signal  Generation  and,  Detection 

It  was  necessary  to  determine  which  method  of 
resonance  excitation  and  signal  detection  was  best  suited 
for  automated  data  acquisition  and  analysis.  The  following 
are  the  three  technigies  s-xaminsd  experimentally: 

a.  swept  frequency  continuous  wave  excitation  with 
tracking  narrowDand  aaterodyne  reception  using 
an  HP  3583A  spectrum  analyzer. 

b.  phase-sensitive  detection  utilizing  a  two-ohase 
lock-in  analyzer  and  a  programmable  frequency 
synthesizer  as  a  signal  source. 

c.  broadband  noise  excitation  in  conjunction  with 
an  FFT  analyzer. 

2-   gat  a  Analysis  Program  Development 

It  was  necsssary  to  determine  a  precise  and 
efficient  computer  algorithm  for  the  fitting  of  the 
resonance  data  to  a  resonance  Lineshape  to  extract  fo,  A, 
and  Q. 

1 1 


3.   System  Implementation 

This  involved  the  assembly  of  the  required 
components  (computer,  analog-to- digital  converters,  counter, 
etc.)  in  a  custom  rack,  interfacing  and  interconnection  of 
components  to  the  computer,  resonator,  and  each  other,  and 
The  writing  of  the  software  :o  illow  the  computer/controller 
to  orchestrate  all  of  the  acquisition,  analysis  and  display 
functions. 

**•   System  2£§£  IskJi  Sli.l!i^ts2S 

The  final  task:  was  to  test  the  system  on  the  lowest 
frequency  plane   wave  modss   of  a   gas  filled,    cylindrical 


resonator 


—  r\  —       —  n      .3  i^  —  <2  — 


etermine:   the  precision  with  which  amplitude. 


Q  and  center  frequency  could  bs  measured;  the  tracking 
capabilities;  the  processing  spesd;  and  the  utility  of  the 
output. 

B.   CONCLUSIONS  FROM  TASK  COMPLETION 

Upon  completing  tasks  1  ani  2,  it  was  determined  that 
method  b  (frequency  synt hesizsr/phase-sensitive  detector) 
provided  the  best  noiss  rejsotion  and  resolution,  was 
computer  controllable,  and  resulted  in  a  data  set  which 
simplified  the  implens nt at  ion  of  the  resonance  fit  algorithm 
chosen  from  task  2.    Method  a  of  task  1  was  rejected  due  to 


12 


the  lack  of  a  suitable  computer  interface  which  would  allow 
automatic  changes  in  frequency  range.  Also,  the  sweep  would 
require  catching  the  values  "on  the  fly"  thus  resulting  in 
inaccuracies  which  were  sweep  rate  dependent  and  data  points 
which   were    not    necessarily    squally    spaced.       Method   c   was    not 


utilized    due  to      the    inability    to    concentrate 


no     - ->  a -r 


nergy   m 


the  bands  cf  interest,  resulting  in  the  ns^d  for  higher 
power  to  the  resor.at.oc  without  an  increase  in  the  signal  to 
noise  ratio  in  the  frequency  donain  of  interest. 

The  completion  of  all  of  these  tasks  resulted  in  the 
system  to  be  describe!  in  the  remainder  of  this  thesis.  The 
system  is  capable  of  a  precision  equal  to  or  better  than  one 
part  in  107  in  frequency,  one  part  in  13,000  in  amplitude, 
and  one  part  in  10  0  0  for  2. 
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II-   ^2MI£^SNT  SETUP  AND  DESCRIPTION 

A.   EQUIPMENT  SETUP  AND  INTERFAZS 

Figure  1  shews  a  block  diagram  of  the  equipment  utilized 
in  the  final  system.  The  following  is  a  brief  look  at  the 
overall  function  of  the  systau  fhich  will  be  followed  by  a 
description  of  each  individual  piece  ci   equipaent. 

The  HP-85  computer  directly  controls  ail  of  the 
equipment,  with  the  exception  :>  f  the  Lock-in  Analyzer,  via 
the  Hewlett  Packard  Interface  Bus  (HPIB}  .  The  computer 
sends  a  value  for  frequency  aid  amplitude  to  the  HP3325A 
Synthesizer/Function  Generator,  which  causes  an  excitation 
of  the  acoustic  resonances  in  the  resonator.  The  amplitude 
analog  signal  receive!  by  the  Lock-in  analyzer  is  converted 
to  a  D-C.  voltage  whish  is  sent  to  the  HP3497A  Data 
Acquisition/Control  tlnit  where  it  is  digitized  prior  to 
being  sent  to  the  zoaputer  via  the  HPIB.  After  processing, 
the  computer  sends  oatput  to  the  HP2673A  Printer  and  the 
HP7470A  Graphics  Plotter.  The  3P3456A  reads  the  temperature 
in  degrees  Celsius  from  a  thermistor  mounted  in  the  body  of 
the  resonator  and,  wnen  interrogated  by  the  HP-85,  sends  the 
current  value  via  the  HPIB.  ,In  the  test  system  (Chapter  IV) 
the   HP3325A  outputs   the  required   frequeacy  and   amplitude 
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signal  by  means  of  a  power  amplifier  to  the  electret 
transducer      at      one      and   of      the      resonator.  An      electret 

transducer  at  the  other  end  of  the  resonator  picks  up  the 
signal  and  sends  it  via  an  amplifier  to  the  E35G  5204.  The 
ZG&G  5204  is  phase- locked  to  tJia  HP3325A  ind  outputs  a  D.C. 
vcitage  which  is  proportional  to  the  Pythagorian  sum  of  the 
amplified  in-phase  aid  quadrature  components  of  the  signal. 
This  analog  amplitude  is  sen-  to  the  HP3497A  which,  when 
interrogated,  sends  via  the  BPIB,  a  digital  measurement,  of 
the  voltage  level.  The  HP-35  also  outputs  values  to  the 
HP2673A    and   HP7U70A    to    be    printed    and    plotted    respectively. 

B.       EQUIPMENT    DESCRIPTION 
1-       HP^SS    Cornea  tec 

The  Hewlett-Packard  Personal  Computer  (HP-85)  is  an 
8  bit  microcomputer  containing  as  standard  a  16K  byte 
memory,  a  built-in  thermal  printer  and  a  small  CRT.  The 
computer      is  programmed      in      3A3I3      computer    language.  As 

utilized   in   this   application,    the    memory   was    expanded   to    32K 
bytes   through      the   use    of   an      add-on    16K   byte      memory   module 
(#52903A).         In   addition,    the    H?IB    interface     (#82937A)  ,       and 
the      ROM      Drawer      (#82  936A)  with      a      plotter/printer      ROM 

(#00085-15002)  and   input/out pat      ROM     (#00085-15003)         were 

utilized. 
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The  HPIB  is  Hewlett-Packard's  implementation  of  the 
IEEE  Standard  488-197  5  and  is  a  parallel  bus  of  16  active 
signal  lines,  3  data  control  lines  and  5  interface 
management  lines.  It  is  oapabls  of  interconnecting  up  to  15 
instruments.  The      input/output    and      prin- er/plott er      ROMs 

further  expanded  the  memory  in  order  to  handle  the 
transformations  necessary  to  transmit  information  via  the 
HPIB   to  the   peripherals. 

2.      aiaj=t£2  Voltaster   HP3455^ 

Although  the  3P3456A  has  numerous  capabilities,  in 
this  application  it  was  only  utilized  as  a  means  -.o  obtain  a 
temperature  from  a  thermistor  and  to  transmit  that  value  via 
the  HPIB  to  the  HP-85  computer.  Its  only  utilization  was  to 
demonstrate  the  system's  ability  to  track  resonances  under 
conditions    of   varying    temperature. 

3-      Sy_n  thesizec^Fi  nction   Generator    HP3  3  25A 

The  HP3325A  has  a  fregusncy  range  from  1  microti  ertz 
to  20  megahertz  vita  a  resolution  of  up  to  eleven  digits. 
The  peak  to  peak  output  amplitude  can  be  set  from  one 
millivolt  to  ten  volts.  It  is  fully  programmable  via  the 
HPIB  and  was  so  utilized  in  this  application.  The  execution 
time    for      a    freguency    command    via      the   bus    is    7.0      msec    plus 
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2.8  msec  for  each  digit  plus  2.9  msec  far  a  decimal  plus 
12.5      msec    for      each      delimiter.  The   execution      time      for 

amplitude    is  6.8    mseo    plus    2.3      msec    for    aach    digit    plus    2.8 
msec   for    a     decimal    piis    13.0    132c    for    a    delimiter. 
u-      Data   Acaui.sition^ContrDL    Unit.    H£3497A 

The  HP3497A  4  as  utilized  as  a  voltmeter  in  this 
application.  It  raoaived  an  aaaiog  voltage  level  fron  the 
lock-in  analyzer  which  was  proportional  to  the  full  scale 
meter  deflection  and  converted  it  into  HPIB  compatible 
signals  to  be  transmitted  to  tas  HP-85.  It  was  completely 
program  controlled,  sampling  and  sending  only  when  so 
directed   by   the    HP- 35. 

5.      SGS.G  Lock-in    Analyzer   ao del   5  20U 

3ecause  tha  5 204  was  not  HPI3  compatible,  desired 
functions  had  to  oa  manually  controlled  vice  program 
controlled.  This  inoluded  changes  to  full-scale  sensitivity 
and      time        constant?.  It      racaived      a        phase      reference 

synchronization  signal  from  tha  HP3325A  and  the  data  signal 
from  the  resonator.  Through  vaotor  manipulation  it  outputs 
the  square  root  of  tha  sum  of  tie  squares  of  the  inphase  and 
quadrature  components  of  tha  amplifiad  signal  to  the 
HP3U97A. 


5.   Power  Am£lifier  HP  a  67  A  and  Signal  Amplifier  HP465A 

These   were   ordinary   amplifiers   inserted   in   the 

system  to   provide  signal  gain.    The   HP457A  was  set   for  a 

voltage  gain  of  two  to  drive   tie  resonator  while  tha  9P465A 

provided  a  gain  of  +23  d3  for  the  output  from  the  resonator. 

7-   In~  el  11  gent  IrajDhics  Printer  HP2673A 

The  HP2673A  prints  bi-iirect ionally  at  the  rate  of 
120  characters  per  s=:or.d  with  a  7  X  1 1  dot  matrix  character 
font  in  a  9  X  15  character  calL  with  twelve  character  sets 
to  choose  from  and  three  modes  Df  formatting.  This  printer 
was  utilized  to  obtain  standard  sized  paper  output  vice  the 
small  sized  thermal  printer  output  of  the  HP-85.  It  was 
utilized  in  the  fully  program  controlled  mode. 

9.   Graphics  Plotter  I?747(H 

The  HP7U73A  was  utilized  to  produce  the  high 
resolution  and  excellent  guality  graphs  included  in  this 
report  and  suitable  for  publication.  It  is  a  two  pen 
plotter  capable  of  utilizing  up  to  ten  colors  through 
programmed  stops.  It  hae  a  ciaracter  plotting  speed  of  up 
to  six  characters  per  second  ia  any  of  five  character  sets 
with  text  written  In  any  direction,  with  or  without 
character  slant,   and  in  varying  sizes.    It  was  utilized  in 
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the  fully  programmed  node  via  t&9  HPIB.  One  of  its  nicest 
features  is  the  ability  to  Safins  the  plotting  area,  thus 
giving  the  operator  tie  capability  of  enlarging  or  shrinking 
the  output  size  as  desired.  AIL  figures  in  this  thesis  were 
produced  on  this  pieos  of  equipnent.  It  is  to  be  noted  that 
although  these  figures  are  actual  system  outputs  or  replicas 
of  system  outputs,  they  were  collected  over  a  period  of  time 
which  involved  several  experiments.  Thus,  each  individual 
variable  may  not  exactly  coincide  in  value  from  figure  to 
figure  since  parameters  suci  as  temperature,  drive 
amplitude,  full  scale  sensitivity,  time  constants,  frequency 
band,  etc.,  ail  affect  eaoh  variable  =nd  were  possibly 
different  for  the  figures  shown. 
9.   Thermistor 

The  thermistor  utiiizsi  in  this  set  up  is  an 
HP-0837-0164  resistive  type  with  a  range  of  -30  degrees 
Celsius  to  150  degrees  Celsius  with  an  accuracy  of  plus  or 
minus  6  degrees  Celsius  for  the  range  of  -75  degrees  Celsius 
to  130  degrees  Celsius.  Beoause  the  temperature  was  only 
being  utilized  as  in  indicator  of  change,  it  was  not 
necessary  to  measure  this  quantity  accurately  for  this 
experiment.      Thus,     the   most    convenient    means    of 
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acccmplishi  r.g  the  iesired  task.  «ias  utilize!.   The  thermistor 
has  a  nominal  resistance  of  53  0  3   ohms  at  25  degrees  Celsius 
rising  to  3.684  Mohms  at  -80   degrees  Celsius  and  falling  to 
92.7  ohms  a t  150  degrees  Celsius. 
10.   Re  senator 

A  detailed  schematic  representation  of  the  resonator 
utilized  to  test  taa  sys-em  is  shown  in  Figure  2  [Ref.  2]. 
Although  this  system  ran  he  a-ilized  with  any  resonator,  a 
description  of  the  one  actually  used  will  be  given. 

The  resonator  is  a  cylindrical  brass  cavity  with 
diameter  and  length  equal  to  2.  54  cm.  The  ends  are  oa oped 
by  simple  eiectret  transducers.  A  small  slot  (0.16cm  X 
1.6cm)  at  the  cylinder  midplaie  exists  to  allow  gas  to  ent=r 
and  leave  the  resonator  (i  heliji  recovery  line)  as  utilized 
in  another  application  [Rsf.  2]. 

A  permanently  polarize!  12  micron  -hick  disk  of 
teflon  (the  eiectret),  aluuinized  on  one  side,  is 
electrostatically  held  against  a  sand  blasted  metal 
electrode,  which  forms  the  electrically  active  element  in 
the  resonator,  while  the  alumlnized  side  is  connected  to 
ground. 
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The  use  of  pr=-charged  alectret  material  eliminates 
the  need  for  bias  supplies  ordinarily  present  in 
conventional  capaciti/e  transdarers.  The  3Dlored-in  circles 
in  the  figure  are  o-rings  utilized  in  assembling  the 
resonator   tc  provide   a    gas    tight    seal. 


23 


III.    ?EQGRA2  EXPLANATION 

A  simplified  £lo«  chart  of  the  program  is  shown  in 
Figure  3.  A  brief  overview  of  the  flow  will  be  given 
followed  by  a  description  of  ths  individual  blocks. 

The  operator  inputs  an  upoar  and  lower  frequency,  a 
drive  amplitude  and  2.  time  constant.  This  •'rime  constant 
will  be  discussed  iacer,  but  i:  also  needs  to  be  set  on  the 
lock-in  analyzer  as  wall  is  sntsred  into  ths  computer.  This 
then  constitutes  tha  input  co  'Search'.  In  'Search1,  an 
approximate  center  irsquency  and  amplitude  is  displayed  for 
all  resonances  within  the  frsguency  band  entered.  The 
desired  modes  are  cidss:.  by  the  operator  and  carried  into 
•Sort*.  In  'Sort'  each  of  the  modes  is  r?measured  and  then 
refined  via  a  five  point  fit  routine  which  results  in  an 
output  of  center  frequency,  maximum  amplitude,  quality 
factor,  and  bandwidth.  A  new  time  constant  is  determined 
and  all  are  carried  into  'Measare'.  In  'Measure',  a  more 
accurate  center  frequency,  maxinum  amplitude  and  Q  are  found 
for  each  mode.  These  then  go  into  'Ravine'  where  they  are 
refined  to  a  high  precision  by  a  modified  Chi-squared 
minimization  algoritha.  The  last  step  is  'Track'  where  each 
of  the  modes  is  continually  checked  as  well  as  the  resonator 
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temperature  until  a  change  in  any  resonant  amplitude  of  10 
percent  or  a  difference  in  temperature  of  0.2  degrees 
Celsius  is  detected  it  rfhich  time  the  program  loops  back 
into  f Measure'  and  repeats  the  routine. 

It   is   appropriate,    at   this   time,    to   state   that 

throughout  the  program,   any  system  equipment  response  times 

(as    mentioned   in    Chapter    II)     necessary   for    data 

transfer/acquisition  are   ailDWrd  to  elapse  through   the  use 

of  WAIT  statements  prior  to  sampling  any  given  response. 

A.   INPUT 

The  operator  selects  a  start  frequency  and  an  end 
frequency  to  define  the  frequency  spectrum  to  be 
investigated.  Also,  a  drive  amplitude  and  maximum  drive 
amplitude  are  selected,  which  can  be  any  value  up  to  and 
including  3500  millivolts  (limit  of  the  HP3325A 
Synthesizer /Function  Generator).  Two  amplitudes  are 
selected  because  later  in  the  program  the  drive  amplitude 
will  automatically  increment  or  decrement  the  value  in  order 
to  keep  the  lock-in  output  in  its  midrange.  The  second 
value  (max.  drive)  is  a  default  value  adopted  should  the 
program  attempt  to  exceed  it  while  incrementing  the  drive 
amplitude.   Finally,  a  time  constant  (tau)  is  selected.   The 
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lock-in  analyzer  divides  the  spectrum  into  256  regions  and 
integrates  over  each.  The  program  calculates  a  maximum  tail 
by 

r=  256/4*  (fuppef  flower) 

in  order  that  there  is  net;  a  spacs  between  regions  and  that 
the  sensitivity  over  the  region  be  nearly  flat.  This  value 
is  a  maximum  such  that  the  operator  selects  a  value 
available  on  the  lock.- in  analyzer  which  is  equal  to  or  less 
than  this  and  enters  the  value  set  on  the  lock-in,  into  the 
computer. 

B.   SEARCH 

The  frequency  range  is  iivided  into  256  parts  as 
mentioned  above.  Each  of  thes=  frequencies  is  sent  via  the 
HPIE  from  the  HP-85  to  the  HP3325A  and  a  resonator  amplitude 
value  for  each  of  these  frequencies  is  obtained  by  the  HP-85 
from  the  HP3497A  via  the  HPIB.  There  is  a  programmed  delay 
(WAIT)  of  4  tau  between  the  sending  of  a  frequency  and  the 
sampling  of  its  amplitude  in  orier  to  allow  the  resonator  to 
attain  at  least  98  percent  of  its  steady  state  value  as 
determined  by  the  look-in  integration  time.  The  choice  of 
256  regions  was  dictated  by  th=  horizontal  resolution  of  the 
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HP-S5*s  CRT  display.  Thase  255  values  ics  then  plotted  on 
the  CRT  with  the  option  of  a  hardcopy  printout  from  the 
HP7470A  as  in  Figure  4.  This  graphically  displays  where 
each  of  the  resonance  modes  within  the  specified  frequency 
interval  is  located. 

From  Figure  4,  it  can  be  seen  that  not  only  are  the 
responses  from  three  major  plaaa  wave  modes  displayed  but 
also  -hose  of  several  other  (azimuthal,  radial,  mixed) 
modes.  It  is  possible  to  conceitrate  on  any  mode  displayed. 
After  the  operator  inters  a  threshold  relative  amplitude 
value,  the  computer  provides  a  printout  on  its  thermal 
printer  listing  all  amplituies  and  their  associated 
frequencies  which  equal  or  exceed  this  threshold  value  (see 
Table  2  in  Chapter  EV  for  printout).  If  the  operator  is 
satisfied  with  the  output,  than  the  plot  and  the  thermal 
printout  are  utilize!  to  input  -co  'Sort1  the  approximate 
center  frequency  and  bandwidth  for  each  desired  mode.  If 
the  operator  is  not  satisfied,  the  program  loops  back,  to 
'Input1,  new  parameters  are  chosen,  and  the  process  repeats 
itself. 
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C.   SORT 

The  purpose  of  'Sort*  is  to  identify  for  the  computer 
the  modes  to  be  investigated  and  to  obtain  "reasonable" 
values  of  center  frequency,   uaxiaum  ampli-ude,   and  quality 

factor  for  each  desired  mods.  "Reasonable" ,  in  this  case, 
is  defined  as  sufficiently  accurate  to  allow  these  obtained 
values  to  be  automatically  controlled  in  'Measure'-  This  is 
achieved  by  the  open  tor  inputting  a  center  frequency  and 
bandwidth  for  each  mode  desired  to  be  worked  upon.  As 
configured,  up  to  nine  modes  can  be  selected.  However,  this 
limit  can  be  increased  through  a  program/dimension  change  to 
whatever  number  is  desired.  Nine  was  chosen  as  a  matter  of 
convenience  to  simplify  memory  nar.agement. 

In  this  portion,  the  bandwidth  of  each  mode  is  divided 
into  one  hundred  (an  arbitrary  number  which  gave  sufficient 
resolution)  equal  frequency  steps  which  are  individually 
sent  and  their  respective  amplitudes  collected  in  the  same 
manner  as  in  search. 

Once  again,  a  plot  is  generated  on  the  CRT  with  the 
option  for  a  hard  copy  (Figure  5)  .  If  the  operator  is  not 
satisfied  with  the  DJtput,  then  the  program  loops  back  to 
the  start   of  'Sort1   and  new  iiputs   are  entered.     If  the 
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operator  is  satisfied,  tha  program  finds  the  maximum 
amplitude  and  its  corresponding  frequency  (center 
frequency)  .  This  is  accomplished  by  examining  the  value  of 
the  amplitude  cf  each  of  the  one  hundred  points  and 
selecting  the  one  which  is  greatest.  The  amplitude  value  of 
the  3  dB  down  points  is  thsn  found  from  this  maximum 
amplitude.  The  program  again  sorts  through  ail  of  the  points 
and  finds  the  closest  "bracketing"  points  to  this  calculated 
3  dB  down  pcint.  k  linear  interpolation  is  then  utilized  to 
calculate  the  corresponding  frequency  for  the  calculated 
3  dE  down  amplitude  point.  Df  course  this  process  occurs 
twice  since  there  is  both  an  loper  frequency  and  a  lower 
frequency  for  the  3  13  down  amplitude  value.  This  process 
constitutes  the  "five  pcint  fit". 

From  these  values,  a  quality  factor  is  calculated  which 
is  equal  to  the  center  frequency  divided  by  the  difference 
cf  the  upper  and  lower  3  dB  iorfi  frequencies.  This  value, 
as  well  as  the  center  frequency  and  maximum  amplitude,  is 
output  to  the  3P-35  thermal  printer  (see  Table  2  of 
Chapter  IV  for  printout). 

This  procedure  is  repeated  for  each  of  the  remaining 
modes.   Figures  6  anl  7  are  the  plotter  outputs  of  the  other 
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two  normal  modes  in ves tiga ted  in  this  example.  Figure  8  is 
an  example  cf  one  possible  reason  •  for  desiring  to  reguest  a 
new  input  (wrong  center  frequeacy) .  In  this  case,  it  is 
possible  that  one  of  the  3  dB  down  points  could  not  be  found 
from  the  data  (cut  of  range) .  From  this  graph,  the  operator 
would  obviously  either  choose  a  new  center  frequency  input 
or  increase  the  bandwidth  suoh  that  the  3  d5  down  point 
would    occur    in    the    data. 


D.        8EW    TAU 

In  this  block.,  three  things  occur.  First,  a  new  time 
constant  is  found.  This  time  iDnstant  is  the  acoustic  time 
constant  reguired  for  each  point  to  attain  63  percent  of  its 
steady      state    value.  This      is    in      analogy      with   the      time 

constant  reguired  foe  a  capaci-or  in  an  electrical  circuit 
to  charge  up  to  its  full  value.  The  acoustic  time  constant 
for  each  mode  is  founl  by  dividing  the  Q  for  the  mode  by  PI 
times  the  center  frequency  of  the  mods.  Because  of  the 
manner  in  which  Q  is  found,  this  eguates  to  dividing  one  by 
PI  times  the  difference  of  the  3  dB  down  upper  and  lower 
frequencies.  After  calculating  the  time  constant  required 
for  each  mode,  the  lew  system  time  constant  selected  is  the 
greatest  one  required,  thus  insuring  -hat  the  worst  case  is 
adequately    covered. 
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Second,  the  background  noise  level  is  determined  by  the 
following  method.  The  synthesizer  outputs  a  frequency  off 
resonance  into  the  resonator  and  the  lock-in  analyzer  is 
sampled     one  hundred      times.  The      mean    and      the      standard 

deviation  of  this  sat  of  measurements  is  calculated.  The 
standard  deviation  is  equate!  to  the  noise  of  the  system. 
It  is  utilized  later  in  determining  the  signal  to  noise 
(SNR)  ratio  for  each  mode.  It  is  also  utilized  in  the 
calculation  of  the  third  item  Df  this  biork  which  will  be 
described   next. 

The  final  (third)  item  is  the  examination  of  each  mode 
to  determine  the  maxiaum  number  Df  points  which  can  be  taken 
without  having  the  possibility  of  the  amplitude  value  of  two 
adjacent  points  being  equal  or  in  the  worst  case  inverting 
position  due  to  associated  uncertainties  or  random  system 
noise.  If  the  frequency  step  size  (  A  f )  is  large  enough, 
then 

AA    =     Amas.    -    Mfo+Af) 

can      be     chosen      so      that      the      inversion      will      not      occur. 
Assuming   a    Eayleigh    listribution    lineshape   then 
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A(f)     = 


A/Q 


V/2 


Z  D 


Replacing   f   with    (fc  ♦»     Af)     ani        6=       Af/fo,      the   resulting 

equation    is 


AA  =  A 


[(2€Q)2+1]1'* 


For  the  system  described,  24  pjinrs  ire  normally  taken  in  a 
frequency  range  of  2  BW  whers  }  =  fo/3W,  Therefore,  the 
following   equation    rasulzs: 

2€Q    =     2  Af/BW    =     1/6    <    1 

and,  that  quantity  sjuarsd  is  <<1.  Expanding  and  dividing 
by    A 


A  A  ,         2  /a/' 
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If   A    is   chosen    to   be    some    number    of    standard    deviations,    &    , 
based   upon    the    noise    measurement,    -hen 


=    *(—QX 


A         SNA 


o 


where      a/k    is    1/SNR    (5 ignal-to-noise      ratio).         Now,       Af   is 

found  tc   be 


"  2(SNl 


fo 


R)  Q 


If   I    is   the    number    of    ncints    -o    bs    taken    in    2    BW    -hen 


2  BW^  2  f0  /8(SNR) 


A  f  Q  A  f         If         N 


If  the  probability  of  an  inversion,  as  described  above,  is 
set  at  0.1  percent,  then  N  =  3.3  based  upon  a  standard 
distribution  curve.  It  is  worth  noting  it  this  point,  -hat 
the  probability  of  ia/ersion  a.wiy  from  maxiaum  is  everywhere 
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lower  until  the  signal  level  is  on  the  order  of  the 
background  noise.  utilizing  the  0.1  percent  criteria, 
I  <  =  1.56*SQR  (SNH>  .  This  number  (I)  is  output  to  the  HP-85 
thermal  printer.  The  program  then  defaults  10  twenty-four 
points  as  a  minimum.  This  value  was  based  upon  a 
calculation  from  the  Nyquist  sampling  rate.  The  minimum 
from  the  Nyquist  sampling  rate  (ir.  integer  value)  would  be 
three;  but  it  is  generally  agreed  upon  in  engineering 
applications  that  a  factor  of  three  or  four  is  required  in 
order  to  obtain  a  satisfactory  reproduction  of  the  signal. 
At  this  point  the  baaiwidth  was  doubled  to  insure  inclusion 
of  the  3  dB  down  points  and  thus  the  ainiraum  number  of 
points  obtained  egualied  24.  For  this  choice  the 
requirement  that  an  Inversion  die  to  noise  nas  a  probability 
of  0.1  percent  or  lees  requires  the  SNE  to  equal  or  exceed 
237  (+47.5  dB)  . 

The  system  has  worked  when  the  indicated  number  of 
points  was  as  low  as  10  (still  lefaulting  to  a  choice  of  24) 
but  occasionally  gave  faulty  values  due  to  the  reasons 
discussed  above.  Although  it  lid  not  always  happen,  it  is 
recommended  that  whenever  the  indicated  number  of  points 
fails  below  24,  eith.ec  the  SNR  is  improved   (increase  drive 
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amplitude,  etc.)  oc  ths  mole  is  net  suitable  fcr  further 
investigation.  A  sample  of  tie  thermal  output  format  for 
tau,  noise  and  number  of  poiits,  is  shown  in  Table  2, 
Chapter  IV. 

2.   MEASURE 

Each  resonance  mole  has  a  frequency  band  calculated  for 
it  in  'Measure'.  This  band  is  twice  the  band  defined  by  the 
3  d3  down  bandwidth.  Also,  the  relative  value  of  the 
amplitude  for  the  center  frequency  is  evaluated  to  ascertain 
the  need  for  adjustment  in  uhs  drive  amplitude.  This  is 
achieved  by  driving  the  resonator  at  the  ranter  frequency  of 
each  mode  and  sampling  the  amplitude  value.  If  it  is 
between  0.3  and  3.95  as  measired  on  the  Data  Acguisiton 
System,  then  no  adjustment  is  necessary.  If  it  is  not 
within  these  limits,  it  is  either  doubled  cr  halved 
depending  upon  whether  it  was  above  or  below  the  limits. 
This  new  value  is  again  tested  in  the  above  manner  until  it 
is  either  within  tie  Units  or  has  exceeded  the  maximum 
drive  amplitude  whioh  was  set  back  in  'Input1.  If  this 
happens  (exceeding)  ,  then  the  system  defaults  to  the  set 
value  of  maximum  irrive  ampli:ade.  This  adjustment  is 
necessary  because  the  full  seals   sensitivity  setting  on  the 
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lock-in   analyzer   is    a    manual   rather   than  computer   controlled 
function. 

In  'Measure*  t ha  resonator  amplitudes  are  obtained  in 
the  same  manner  (frequency  sent--amplitude  measured)  as  in 
•Sort',  only  this  time,  twanty-four  freguencies  vies  one 
hundred  are  utilized  and  tha  'WAIT'  becomes  12  tau  which 
allows  the  system  :o  reach  within  6  ppm  of  the  steady  stats 
value  to  fully  utilize  the  5  1/2  digit  A-to-D  conversion  of 
the  H?  4397A  (as  explained  previously  at  the  end  of  'NEW 
TAUf).  Once  the  data  are  collected,  the  following  seven 
point  fit  is  applied.  The  maximum  amplitude  point  and  the 
closest  adjacent  points  on  either  side  of  it  are  used  to 
define  a  parabolic  :arve.  Iha  maximum  amplitude  of  this 
parabola  and  the  corresponding  center  frequency  are  found  by 
setting  the  first  derivative  of  the  amplitude  with  respect 
to  frequency  to  zero.  Solution  of  the  resulting  equations 
and  subsequent  insertion  of  tha  solution  into  the  original 
equation  yields  a  unique  laximum  amplitude  and  center 
frequency.  This        amplitude      is      then      utilized        in      the 

calculation      of      the      3    dB      dowi      amplitudes.  Frequencies 

corresponding      to  tha       3dB      down    points      are      then    found      by 
linear      interpolation      in      the      same      manner      as      previously 
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described  in  'Sort'.  A  new  juality  factor  is  calculated 
utilizing  these  new  frejuencias  and  all  three  of  these 
values;      2,      maximum  i mplituda    ind    center    frequency   are    sent 

into    'Ravine*. 

F.        RAVINS 

The  values  for  maximum  amplitude,  center  frequency  and 
quality  factor  determined  by  the  seven  point  fit  described 
in  the  previous  section  ('Measure')  are  substantially  more 
precise  than  those  possible  using  conventional  techniques 
employing  a  wave  analyzer  or  spectrum  analyzer  under  manual 
control.  However,  the  high  speed  numerical  processing 
ability  of  the  digital  computer  allows  a  greatly  expanded 
treatment  of  the  data  for  Daly  a  minimal  investment  in 
actual  processing  tims. 

Implicit  in  ths  assignaeit  of  A  (amplitude) ,  fo 
(freguency),  and  Q  (quality  factor)  as  the  three  parameters 
which  characterize  a  resonance  line  shape  (as  is  done 
throughout  this  thesis)  is  the  assumption  that  the 
amplitude,  as  a  function  of  freguency,  near  resonance  can  be 
described  by: 
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A(fn)  =  A/Q  [(fn/fo-fD/ea)  2  ♦  (Vg)2]1^ 

where  :   A  ( fn)  is  tha  r?£DP.an:5  lineshape 

A  is  the  calculated  maximum  amplitude 
Q  is  the  calculated  quality  factor 
fo  is  the  calculated  canter  frequency 
fn  is  the  maasurad  fra:juancy  (24  par  mode) 

Using  all  twenty-four  points,  it  is  possible  to  obtain  the 
sum  of  the  squares  of  deviations  of  the  data  from  the 
resonance  lineshape  dafirad  by  a  particular  choice  of  A,  fo 
and  Q  as: 

24 


S2  (A,fo,Q)  =  ]^An  "  *  (fn)] 


n=l 

where    An    is   the    measured   amplitude      of   the    n  (th)       data    point 
and   A(fn)     is  as    abova. 

Conventional  data  analysis  then  accepts  the  "best 
choice"  for  A,  Q,  and  fo  as  that  which  causes  S2  to  be  a 
minimum  [Ref.  3  ].  This  minimization  for  a  linear  function 
is  a  straight  forward  algebraic  process.  However,  an 
analytic  least  squars  fit  to  a  non-lineac  function  requires 
matrix  inversions.  This,  in  turn,  not  only  requiras  large 
blocks  of  memory  and  computer  time,  but  also  can  easily  lead 
to  errors  due  to  internal  round-cff  or  truncation  of  digits 
within      the      computar    [Ref.    4].         For      these      reasons,         the 
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following  method  was  devised  and  utilized  for  each  of  the 
resonant  modes  being  i nvestigated.  Briefly,  S2  was  computed 
for  small  variations  in  A,  fo  and  Q  as  well  as  for  the 
originally  computed  values  for  these  terns  from  the  seven 
point  fit.  These  S2  terms  were  then  utilized  in  a  parabolic 
fi-1:  in  an  attempt  to  find  those  values  of  i,  fo,  ar.d  Q  which 
wculd  cause  S2  to  be  a  minimua.  Because  the  results  of  the 
seven  point  fit  were  already  quite  good,  the  procedure, 
briefly  described  above  and  expanded  upon  in  the  remainder 
of  this  section,  ooiverges  ver/  rapidly  to  an  astonishingly 
high  degree  of  precision.  This  precision  will  be  analyzed 
and    discussed    in   detail    in    Chap*:  si:    V. 

For  the  remainder  of  this  paper,  the  sum  of  the  squares 
cf  deviations  (S2)  fill  be  referred  to  as  the  C-G  value, 
which   was    chosen   for    historical    reasons. 

One  at  a  time,  each  of  the  three  parameters  (A,  fo,  Q) 
in      the      resonance      LLneshape      (A(fn))         are      varied.  The 

following   is  an   example   to   clarify   this   process. 

First  the  C-G  tern  for  the  calculated  2  is  found.  Mow  Q 
is  varied  by  +0.5  peccent  and  a  second  C-G  term  is  found. 
Finally,  2  is  varied  by  -3.5  percent  and  a  third  C-G  term  is 
found.      These    three    C-G   terms    ace    used   to    define    a    parabolic 
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curve  (Figure  9)  and  a  new  value  for  Q  which  minimizes  the 
C-G  term  is  found. 

This  new  value  3f  2  replaces  that  found  in  the  seven 
point  fir  and  becomes  tne  new  value  to  be  inserted  in  the 
resonance  iineshape.  This  tine,  however,  A  is  the  variable 
to  be  adjusted.  Just  as  occurred  with  2,  C-G  terms  are 
found  for  A,  A  +  3.27»h,  and  A  -  0.2%A.  Again,  these  C-G 
terms  are  utilized  in  a  paraaoLic  curve  minimization  and  a 
new  value  for  A  is  obtained  which  now  replaces  the  A  found 
in  the  seven  point,  fit.  Finally  -hese  new  values  of  Q  and  A 
are  utilized  to  determine  a  new  fo  in  the  same  manner.  This 
completes  one  iteration  of  the  process  and  thus  provides 
these  new  values  3f  2,  fo,  aid  A  for  the  next  iteration. 
Each  iteration  refines  the  precision  of  the  values  of  Q,  fo, 
and  A  and  provides  this  new  set  of  values  for  the  next 
iteration.  This  process  is  repeated  until  the  next 
iteration  contributes  a  negligibly  small  variation  in  the 
C-G  value  [  Hef .  3]. 

In  the  actual  program  utilized,  the  terms  were  variei  in 
exactly  the  order  described  because  in  the  relative 
precision  of  the  terns,  Q  was  tie  least  precisely  known  term 
(historically,   this  is   always  true).    Since  this   was  the 
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case,  varying  it  first  resulted  in  the  greatest  effect.  The 
next  icwest  precision  term  wis  A,  so  it  was  obviously  the 
next  choice  to  be  varied.  Finally,  fo,  which  was  actually 
well  known,  was  varied.  Beoaise  of  the  relative  precision 
in  Q,  A,  and  fo  the  variatiois  were  arbitrarily  chosen  at 
the  indicated  values  (0.53  zf.  2,  0.21  of  A,  and  0.035^  of 
3»)  ,  as  a  matter  of  experience  siivce  these  values  reflected 
a  sufficient  variatian  for  the  'RAVINE1  technique  to  be 
efftctive.  The  better  the  approximations  for  Q,  A,  end  fo 
put  into  'RAVINE*,  the  smaller  these  values  car.  be  made. 

Upon  completion  of  the  final  iteration  (program 
dependent)  the  three  newest  values  for  3,  fo,  and  A  are 
utilized  to  find  a  final  C-G  term  which  is  then  printed  on 
the  system  printer  aid  stored  31  a  data  tape.  This  C-3  term 
is  related  to  Chi-sjuared  ani  may  therefore  be  usei  to 
determine  the  "goodnsss  of  fit"  [  Ref .  3  ]•  At  this  time, 
other  parameters  and  variables  are  also  printed  on  the 
system  printer  and  stored  on  the  data  tape  (time, 
temperature,  pressure,  moie  number,  fo,  A,  Q,  SNR,  drive 
amplitude,  sequential  data  tape  item  number,  and  C-G  Ravine 
value).  These  will  be  discussed  at  the  end  of  Chapter  IV 
and  a  sample  output  will  be  shown  (Table  1)  . 
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If  there  are  additional  aoies  to  be  evaluated,  then  a 
loop  is  made  bask  to  ' Measure1  and  the  next  mode  is 
evaluated  in  the  sane  manner  as  described  above.  Upon 
completion  cf  the   final  mods  to  be   evaluated,   the  program 

shifts  into  'Track'. 

G.   T3ACK 

In  this  portion  the  currently  measured  amplitude  of  each 
mode  is  compared  with  the  isasured  value  obtained  during 
'Measure*.  If  the  value  differs  by  +/-  10  percent 
(indicating  a  change  in  some  external  parameter),  then  the 
program  loops  back  to  'Measure'  to  find  the  new  Q,  A,  and  fo 
for  each  mode.  Also,  the  temperature  of  the  resonator  is 
measured  and  if  a  difference  of  0.20  degrees  Celsius  is 
found  then  the  system  returns  to  Measure'.  If  neither  the 
amplitude  nor  the  temperature  has  met  the  criteria  necessary 
to  return  the  system  to  'Measure',  then  the  'Track'  section 
continually  loops  through  itself  until  either  cf  these 
criteria  is  met.  The  system  remains  in  this  endless  loop  of 
'Measure',  'Ravine',  and  'Ira-k'  until  the  operator  causes 
the  system  to  halt,  the  lata  overflows  the  available 
storage  space  on  the  data  tape  (850  node  entries)  or  a 
catastrophic  system   failure  occurs  (loss  of   power,   etc.). 
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The  preferable  methyl  of  stoppage  is  to  utilize  a  'PAUSE1 
command  from  the  computer  Icsyboarl  siice  this  enables 
continuation  or  program  restart  as  desired.  It  is  to  be 
nczed  that  any  keyboard  action  will  result  in  a  halt. 
This,  as  well  as  any  of  the  other  mentioned  halts  (except 
'PAUSE')  may/will  regiire  re-loading  the  program  into  memory 
prior   to   re-running    the   system. 
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17.        SAMPLE    EXPERIMENT 

In  order  tc  test  the  system's  ability  tc  measure  the 
stated  quantities  and  to  track  the  resonant  modes  as  some 
external  parameter  was  varied,  a?,  air  filled  resonator  was 
cooled  in  a  liquid  nitrogen  bath.  This  system  was  chosen 
for  Its  simplicity  aid  because  the  variation  in  sound  speed 
with      temperature   foe      an      ideal    gas      is      well    known.  The 

resonator  was  placed  in  a  Dewar  flask  and  a  small  quantity 
of  liquid  nitrogen  <is  introduced  to  cool  the  resonator. 
Once  the  temperaturB  staDilizei,  the  computer  system  was 
activated  and  the  lowest  threr  plane-wave  resonance  modes 
were  tracked  as  tai  system  warmed.  This  tracking  and 
recording  of  data  proceeded  unassisted  by  the  operator  for 
approximately  twenty-two  hours.  At  the  end  of  this  period, 
the  resonator  had  reached  room  temperature  and  the  process 
was  halted.  The  printout  was  examined  and  the  data  tape 
was  utilized  to  produce  Figures  10  through  14  en  the  system 
plotter   to    demonstrate    the    utility   of    the   system. 

Figure  10  shows  a  point  by  point  plot  of  temperature  in 
degrees  Kelvin  versus  time  in  seconds.  After  time  10,000, 
the  curve  begins  to  appear  as  a  solid  line  because  the 
individual    points    become   too    closely      grouped    together    to    be 
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visually  distinguished  in  the  pLot.  This  likewise  occurs  at 
varying  places  in  the  other  figures.  Also,  in  Figure  10r 
there  is  a  slight  dip  at  the  beginning  of  the  curve  because 
the  system  was  activated  as  the  temperature  was  still 
slightly  decreasing  vice  stable  as  descriDed  aocve. 

Figure  11  shows  a  plot  Df  the  log  of  the  difference 
between  a  reference  temperature  (a:  infinite  time)  and  each 
individually  measured  temperature,  versus  time.  This  plots 
as  a  nearly  straight  line  whoss  slope  is  a  measure  of  the 
exponential  thermal  sjuilibriuia  time. 

Figures  12,  13,  and  14  show  plots  of  the  square  of  the 
normalized   frequencies  of   eacn   mode   (fo/M#)  2  versus   the 

aas  these  olots  should 


absolute  terrroeratur 


ror  an  loaa 


be  straight  lines.  In  a  perfszt  resonator,  each  successive 
resonant  mode  would  oe  an  exaot  multiple  of  the  first  and 
the  normalized  freqisncies  would  lie  identically  upon  one 
another.  In  this  experiment  node  two  is  parallel  to  mode 
one  but  slightly  higier  in  normalized  value  while  mode  three 
lies  diagonally  in-between  mode  one  and  mode  two.  Because 
the  exact  physical  properties  of  the  resonator  itself  were 
not  being  investigated,  this  was  not  persued  further.  The 
slight  curvature   of  the  lines  is   probably  due  to   the  fact 
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that,  the  thermistor  ioes  not  indicate  absolute  temperature 
exactly.  The  slop?  of  the  Lines  and  the  length  of  the 
resonator  could  be  utilized  to  determine  the  constants  in 
the  ideal  gas  equati^i  of  state. 

These  figures  (13  through  "U I  demonstrate  a  few  of  the 
possible  means  of  displaying  the  data  stored  on  the  tape. 
It  is  intended  that  future  users  will  determine  what  plots 
are  needed  and  generate  them  accordingly.  The  program  used 
to  produce  one  of  these  plots  irsrn  the  data  tape  is  included 
in  Appendix  3. 

To  enumerate  what  data  is  available  on  the  tape,  a 
partial  reproduction  (dot  matrix  output  is  not  permitted  in 
this  thesis)  of  the  system  printer  output  for  another 
experiment  is  shown  in  Table  1.  The  time  shown  is  the 
average  time  at  whioh  the  measurements  were  taken.  The 
temperature  is  the  average  at  which  time  the  measurements 
were  taken  (in  degrees  Celsiusi  .  The  pressure  is  set  to 
zero  because  a  varyiag  pressure  was  not  to  be  evaluated,  and 
thus,  a  means  of  detecting  pressure  was  not  included  in  the 
equipment  setup.  The  next  entry  in  the  line  is  the  mode 
number  for  which  the  data  applies,  followed  by  the  final 
ravined  values  for  fo.  A,   ani  2-    The  SNR  is  obtained  from 
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the   division  of  A    by   the  value   found  for   the  noise   (as 

explained  in  Chapter  III)  .  The  next  value  on  a  normal 
printout  would  be  the  drive  ampLitude  of  the  system  followed 
by  the  sequentially  increasing  umber  of  the  data  tape  item. 
These  two  items  were  deleted  from  this  figure  in  orler  to 
conform  to  thesis  requirements  without  requiring  photo 
reduction.  The  whole  line  constitutes  one  lata  item  so  far 
as  the  -ape  limit,  of  350  lata  itsms  are  concerned.  However, 
each  item  in  the  line  is  an  actual  separate  entry  within 
each  of  the  data  it  en  5.  The  final  entry  is  the  C-G  'Ravine* 
value. 

Table  2  is  a  reproduction  of  the  output  obtained  from 
the  HP-85  thermal  printer  iuring  ths  course  of  an 
experiment.  For  farther  amplification  of  each  of  the 
entries,  the  reader  is  directed  to  Chapter  III,  sections  B, 
C,  and  D.  The  system  printer  was  utilize!  for  this  figure 
again  because  dot  matrix  outputs  are  not  permitted  in  this 
thesis. 


50 


TABLE  2  SAMPLE  OUTPUT  OF  HP-85  THERMAL  PRINTER 


AMPLITUDE  IN  VOLTS   FREQUENCY 


01082675 
00829425 
01056175 
00820725 
00828725 
01466500 
02054900 
01539075 
00869400 


6817.00 
6896.  00 
13769.  00 
13848.  00 
20405.  00 
20484.  00 
20563.00 
20642.00 
20721.  00 


MODE  1  CENTER  FREQ  IS  6830  AND  AMP  IS  .0181355  Q  IS  19a 123117729 
MODE  2  CENTER  FREQ  IS  13756  AND  AMP  IS  .0127195  Q  IS  135.493980736 
MODE  3  CENTER  FREQ  IS  20540  AND  AMP  IS  .02531475  Q  IS  194.224445807 

TIME  CONST  >  =  a  86062488511   NEW  TIME  CONST  IS  10 

The  mean  is  .000378561   The  standard  dev  is  9.  79752854S4E-7 

NO.  OF  POINTS  FOR  MODE  1  IS  212  NO.  CHOSEN  IS  24 
NO.  OF  POINTS  FOR  MODE  2  IS  177  NO.  CHOSEN  IS  24 
NO.  OF  POINTS  FOR  MODE  3  IS  250     NO.  CHOSEN  IS  24 
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?•   EVALUATION  AND  SECOMMENDAI?  IONS 

This  section  describes  the  method  utilized  to  determine 
the  "goodness  of  fit"  and  the  -resulting  precision.  The 
program  previously  described  was  utilized  in  a  steady  state 
condition  (no  variation  of  estsrnal  parameters).  Mso, 
•Ravine'  was  programmed  for  fifteen  iterations  per  mode, 
with  a  printout  of  the  0,  A,  and  fo  values  for  each- 
iteration.  In  this  nanner  the  change  per  iteration  could  be 
evaluated  for  each  of  the  variaolss.  As  seated  previously, 
'Ravine*  is  completed  (Z-G  minimized)  when  successive 
iterations  produce  a  negligible  variation  in  the  value 
obtained  [  Hef .  3  ]•  ?  or  a  given  sample  run,  the  variation 
between  the  fourteenth  and  fifteenth  iteration  ((15-14) /15) 
was  approxi  nately  ons  in  one  miLlion  (mode  1)  to  one  in  ten 
million  (mode  2)  for  2,  one  in  two  million  (mode  1)  to  one 
in  fifteen  million  (mode  3)  for  A,  and  one  in  six  billion 
(mode  1)  to  one  in  one  hundred  thiry-seven  billion  (mode  2)  . 
It  is  evident  that  tie  change  per  iteration  is  definitely 
small  and  thus  the  fifteenth  iteration  values  were  utilized 
as  a  base  line  for  tia  calculation  of  precision  to  follow. 
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TABLE  3     COMPARISON  OF  PRECISION  OF  ITERATIONS 


MODE  I 

MODE  2 

MODE  3 

Q 

Q 

Q 

3rd/ 15th 

1.  04E-3 

8. 94E-5 

1.  25E-4 

5th/ t 5th 

3.  85E-4 

3.  31E-5 

4.  64E-5 

3rd/ 1 5th        4. 66E-4      2. 35E-5     3. 99E-5 
5th/15th        1. 72E-4      1.C5E-5     1. 47E-5 


3rd/15th        1. 20E-7      3. 62E-9     1. 54E-8 
5th/15th        4.45E-8      1. 34E-9     5. 71E-9 
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Each  C-G  ravine  value  for  the  third  and  fifth  iteration 
was  compared  to  the  base  line  to  determine  the  precision  of 
the   system.        This      mis   accomplished    by   taking      the   ratio   of 

the   difference      of   tha      base   Una    and      the    iteration      to    the 
base    line    value.  Table    3    shows    the      approximate    (rounded) 


values      obtained     for      each   by      node.  This 


i.  tMicoeuii 


best/worst  precision  of  one  in  ten  thousand/one  in  one 
thousand  for  Q,  of  one  in  thirty  thousand/one  in  two 
thousand  for  i,  and  one  in  thraa  hundred  aillion/one  in  one 
million  for  fo  for  the  third  iteration.  The  program,  as 
presented,  stops  aftac  the  third  iteration  because  each  full 
iteration  requires  45  seconds  on  the  HP-85  and  that  was  the 
point  chosen  in  tha  precision/time  trade-off  for  this 
particular  application.  This  will  be  further  discussed  at 
the  end  of  this  chapter.  The  results  for  the  fifth 
iteration  (in  the  same  manner)  are  one  in  thirty 
thcusand/one  in  three  thousand  for  Q,  one  in  one  hundred 
thousand/one  in  five  thousand  for  A,  and  one  in  five  hundred 
million/  one  in  twenty  million  for  fo. 

Precision  has  bean  discussei  since   it  is  independent  of 
the   equipment  utilized.     Accuracy,    on   the  other   hand. 
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involves  net  only  the  precision  to  which  any  given 
instrument  measures  as  well  as  -hat  of  the  Dverall  system 
but  also  the  accuracy  of  each  aeasuremer.t  and  the  accuracy 
of  the  overall  system.  For  example,  this  system  will 
measure  the  given  value  of  the  temperature  to  within  mere 
than  five  significant  figures  between  -90  degrees  and  150 
degrees  Celsius,  yet ,  it  is  3nly  accurate  to  within  six 
degrees  Celsius  (a  fuction  of  the  thermistor).  Thus,  there 
is    a    dramatic    difference   between    precision    and    accuracy. 

A.       CONCLUSIONS    AND    RE COMHSNDATIONS 

The  ability  of  this  systetn  to  obtain  2,  A,  and  fo  far 
surpasses  that  of  any  manual  method.  In  aftiition,  it  allows 
the  opera- or  to  perform  additional  functions  as  the 
experiment  progresses,  dees  not  make  transcription/ 
interpretation  errors  as  is  pDssible/liicely  with  manual 
recording  of  data,  gives  highly  reproducible  results 
independent,  of  the  skill  of  the  operator,  and  allows  for  the 
resulting  data  to  oe  to  be  easily  presented  in  a  variety  of 
formats   in    a  very   short   perioi    :>£    time. 

It  was  originally  intended  that  this  complete  system  be 
a  prototype  to  then  b2  transferred  to  a  system  utilizing  the 
HP9836      rather    than      t he      HP-35.  However,       the      equipment 
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arrived  too  late  to  be  incorporated.  One  major  advantage  of 
the  HP9836  is  that,  it  has  an  5ight-to-or.e  speed  advantage 
over  the  HP-35.  Speed  is  the  reason  that  only  three  rather 
than  five  iterations  of  'Ravine'  were  utilized.  Sines  each 
mode  took.  apprcxiaat  sly  135  seconds,  6.75  minutes  elapsed 
from  the  start  cf  *  Measure1,  Sole  1 ,  to  the  next  start  of 
'Measure',  Mode  1.  Obviously,  tha  greater  the  number  of 
modes  the  greater  tha  time  required.  This  time  diffferer.ee 
could  be  a  major  problem  if  tha  changing  external  parameter 
caused  the  resonance  of  the  uoles  tc  shift  significantly  so 
that  the  response  no  longer  fitted  within  the  chosen 
bandwidth.  One  method  of  overcoming  this  would  be  to 
incorporate  a  prediction  of  the  next  location  based  upon  the 
change  of  the  parameter.  This  was  not  incorporated  in  this 
program  for  two  reasons.  One,  it  was  never  necessary  in  the 
experiment  run;  and  two,  the  shift  to  the  HP9836  with  its 
speed  advantage  wo  till  eliminate  this  lengthy  time  period  (it 
could  perform  twenty-four  aoles  in  the  same  period  it 
presently  takes  for  three)  . 

Another  future  change  couli  be  to  work  in  'Ravine*  for 
the  first  set  of  data  while  taking  the  measurements  for  the 
second  set  cf  data.    This  wouii   be  a  better  utilization  of 
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computer  tiros  since  t hers  are  currently  periods  when  the 
computer   is    idle    (tfAir    periods    and    time   constants). 

Another  improvement  would  be  to  replace  the  lock-in 
analyzer  with  a  computer  controllable  model.  This  is,  of 
course,  a  trade-off  of  price  versus  convenience/time  sines 
there    is   approximately    a   two-to-one   ratio   in    equipment   cost. 

There  are  also  nunercus  coaputer  programming  techniques 
which  could  shortsn  portions  of  ths  program  and/or 
associated  run  times.  One  obvious  change  would  be  to  make 
CRT  presentations  of  graphs  optional  or,  if  the  system  under 
test  was  already  well  known,  to  -nake  the  option  available  to 
enter  the  program  at  the  'NEW  r&0*  point  thus  eliminating 
the  need  to  establish  the  desirai  mode  data.  This  was  not 
done  in  tne  presented  program  since  it  will  be  utilized  by  a 
variety  of  users  and  a ach  will  be  able/is  encouraged  to  make 
modifications    for  spsoifio    applications. 


APPENDIX    A 


REFERENCE    GUIDE    TO    PROGRAM    LISTING 


Line    Nos. 


D=  scription 


1    -    251 

270    -    640 
650    -    1690 
1691    -    1899 
1900    -    2676 
2679    -    U997 
4998    -    5320 


r  s  ?  ut 

5  E  AP.  CH 
SORT 
NEW    TAU 
aEASORE 

RAVINE 
TRACK 


5  8 


1  DIM    F4(200)  ,A <3G0) ,C( 100) 

2  CLEAP 

3  PRINTER    IS    2 

4  OuTrLT    7  22    ; "F4M6" 

5  DISP    "REMOVE    PROGRAM   TAPE    AND    INSTALL    DATA    TAPE"     fc    EEE 
P 

6  LISP    "v\riEN    COMPLETED  ,    ENTEP    1" 

7  INPUT    Y 

6  IF  Y=l  THEN  S  ELSE  CLEAR  G  EEEP  @  CC1C  2 

9  ERASETAPE 

10  CREATE  "DATA"  ,850  ,88  @  REWIND 

11  CLEAR  e?  EEEP  i  intlai  sec  up 

13  SETTIME  U  ,0 

2  1   LISP    "Enter    lcwer    rreq" 
31    INPUT    Fl    i    lcwer    free 
41   CLE^R    fe    EEEP 
51    LISP    "Enter    upper    rrec" 

01  INPUT    F2    !    upper    rrec 

7  1  IF  F2<  =  ri  THEN  CLEAR  @  BEEP  £  GOTO  51 
61  1=1/ ( 4* ( (F2-F1) /25b) )  !  time  constant 
SI    CLEAR    @    BEEP 

101   DISP    "Set    time   constant   on    5204    at    ";T;"    or    smaller" 

111   DISP    "When   complete,    enter    value    set    (in    ir  illi-sec)  " 

121    INPUT    Tl 

131  IF  T1>T*1000  THEN  CLEAR  g  BEEP  @  GOTO  101 

14  1  CLEAR  H    EEEP 

151  DISP  "Enter  amp  cr  criv.ing  rrec  in  rr.v  (.=<3  50G)RMS" 

161  DISP  "DECIMAL  VALUES  ARE  NOT  PERMITTED" 

171  INPUT  A  !  airp  of  rreq 

131  IF  A>3500  THEN  CLEAR  @  EEEP  @  GOTO  151 

18  2  FOR  N»l  TC  S 

18  3  A8(N)*A 

18  4  NEXT  N 

191  CLEAR  (§  EEEP 

201  DISP  "Enter  largest  value  cr  amp  ever  oesireo  (=<350 

Gmv  RMS)  " 

211    INPUT   Al    !    ii.ax    future    amp 

221    IF    Al>3500    THEN    CLEAR    @    BEEP    %    GOTG    201 

2  30    !    these    are    tne    values    tc    set   up    tr.e    3325A 
231  A$=vAL$ (A) 

2  32    A2$»"AM" 

23  3   A  3 5= "MR" 

240    OUTPUT    717    ;A2$,A$,A3$ 

24  1    r.\AIT    20  0 

251  F  3=CEIL  i  (F2-Fl)/256)     i    smallest    int    >=   bandwidtn 

252  !  SEARCH 

25  3  !  SEARCH 
254  !  SEARCH 

69 


it 


33G 

350 


1 


■^  -j  £. 

3o  0 


261  i  next  grcup  sencs  free,  tc 

34  9  7A 

26  9  CLEAR 

2  7 0  DISP  "I  arc  working  in  SEARCH" 

271  FOR  N»0  TO  256 

2  51    F4 (N) =Fl+N*F 3 
2o2    1    F  4  ( N )     is    rrec    sent 

290  F4$=VAL$  (F4 (N)  ) 

291  H1?="HZ" 

292  H2$="FR" 

3 Cu    OUTPUT    717     ;K2$fF4$,tfl$ 
31b    wAIT    4*T1+10C 

OUTPUT  7C9  ;MvT3" 

ENTER  70S  ;  A  ( N ) 

NEXT  N 

!  F5  ANL  tb    ARE  SCALE  FACTORS 

EEEP 

F  5=f  1-. 1* (F2-F  i) 
370  F  b=F2-H  .  1*  ( E  2,-F  I ) 
360    COLLAR    (s    CLEAR 

3  SO    SCALE    F  5  ,F  c  ,-  .  1  , 1  .  2 
4 CO    XAXIS    Or1000#Fl,F2 
4 10    JfAXIS    Fl  ,  .  1  ,-.1  ,1.2 

4  11    FOF.   N=0    TO    2  5b 
4 12    REN OR 

4  13    PLOT    F4(N)  ,A (N) 
4  14    NEXT    N 

4 1 5  C I S  P    "IF    Y  0  b    D£S I RE    - 
SEEP 

416  INPUT    Y 

41  "3    IF    Y=2    THEN    47  7 

4 10    CCLEAR    §    CLEAR    £    EEEP 

4  19    CISP    "ENTER    FIGURE    #    TO    cE    PRINTED" 

420  INPUT    L$ 

421  PLOTTER    IS    7C5 
4  22    PEN    1 

423  SCALE    t  1-  .1*  (F2-F1)  ,F2+50G  ,-..3  ,1  .  3 

424  .XAXIS    0,1000  ,Fl,F2 

42  5  YAaIS  Fl  ,. 1,0  ,1.2 
4  2  6  PEN  2 

42  7  FOR  N=u  TC  2  5b 

42  6  PEN OP 

42  9  PLOT  F4(N)  ,A  (N) 

4  30  NEXT  N 

431  PL NOP  iy  PEN  1 

4  32  LOIR  0, SIN  (SO.)  @  PENOP 

433  FOR  X=F 1  TO  F2  STEP  2  0GO 

4  34  MO  vE  X  ,- . 1 3 

4  35  LALlL    vALb(X) 

436    NEXT    X    @    PENOP 


5A  ,    ana   qets    arc: 


rcrr, 


Ct'i     EN  TER 


EL 


C.  Ja  J.  £  R 


(5    PtNOP 


70 


440  LDIR    0 

441  FGE    Y=o    TC    1.2    STEP     .1    fe    PENUP 
4<*2    MOVE    F  1-.  05*  (F2-F1)  ,Y 

44  3    LABEL    VAL$(Y) 

44  4  NEXT  Y  @  PENUP 
4  5u  EC  IE.  C  g  PEN  LP 

4  51  MO v £  F1+30G0,-.21 

452    LABEL    "FRECjUENC¥    IN    KZ" 

4  52    LCIE    0    is    PEN  LP 

45  4    i-'C  vE    F  i  ,-  .  2 

4  55    L^tuL    "PLGl    Cr     RELATIVE    AMPLITUDE    VS    FE.EC;U£NCY     (FULL 

SFECTROM  " 
45b    LLIR    0, SIN (90)     't    PENUP 
457    NC  v  £    FL-.uc*  |r2-tl)  ,  . ': 
Oo    LABEL    "RELATIVE    AMPLITUDE" 
4  5  9    LDII     L    is    PEN  LP 
^b G   MOVE    F I +2 0  00  ,1 .2 
461    LALEL    "FIGURE    ";L$ 
4b2    PENUP 

4  7b    PLOTTER    IS     1 
47  7    CLEAP     is    o££P 

47S    DISP    "Lc    ^cu   cesire    tc    rerun    exp    with    cirierer.t    ura 

ir.eters" 

47 9  CISP  "I i  sc  enter  1,  ctnerwise  enter  2" 

^£0  INPUT  Y 

490  GCLEAR  ^    CLEAR  @  BEEP 

5  0  0  IF  Y  =  1  TH  £  N  1 1 

51C  DISF  "Enter  cecisicn  point  rcr  cirpiituce  (0  TC  1.2)" 

520  INPU1  A4  i  CECISICN  POINT 

5:0  CLEAR  e  EEEP 

54o   CISP    "Enter    rull-scale    sensitivity    setting    from    tne 

52u4    in    volts" 

550  INPUT  A5  !  SENSITIVITY  SETTING 

56  0  PRINT  "Amplituce  in  vclts" ; "    Frequency" 

5  70    PRINT 

580    FOE    N=0    TO    2  5o 

590    IF    A(N)<A4    TEEN    6  30 

6  0b    A  (N) =A5*A(N) 

610    PRINT    USING    620    ;    A(N)  ,F4(N) 

o20    IiVjAGE    IX  ,  L. CCLECCCC , 10X ,CDCLCC .LD 

6  30    NEXT    N 

631  PRINT  USING  6  22 

o22  INAGE  3/ 

640  CLEAR    @    £>EEP 

641  !  SORT 
6  42  !  SORT 
64  3    i    SORT 

6  50   CISP    "hew   many    ir.cces    <uc  ycu    cesire    tc    track    (i-AX    CF 
9)" 
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o  o  6 

c7o 

660 

o96 

700 

710 

720 

7  30 

740 

750 

7o0 

761 

ng 

770 

7  7  L 

7c  0 

"So 

SCO 

610 

6  20 

o  30 

631 

c40 

656 

651 

co  o 

670 

obO 

6  90 

i?66 

Sic 

S2G 

930 

*4G 

941 

^42 

EF 

94  3 

*44 

9*5 

946 

947 

954 

y55 

956 

957 

956 

959 

96u 

INPUT    fc!     !    NOMBEP    OF    MODES 
r  0  F    N  =  1    i  O    M 
CLEAR    t    sEEP 

DISP    "Wfcat    is   center    rreq    tor    incde    ";N 
INPUT   fc(N)     !    CENTER    FEEC 
CLEAR    @    EEEP 

DISF    "vvhat    is    rrecwjotr.    rcr    n.cce    "  ;K 
DISP    "MUST    tt    CR EATER    THAN 
INPUT    C(N)      !     FREcj-tilDTH 
N  E  X  i     In 
CLEAR 

!    Frcn,    77u    tc    ^oO    is    the    first    rcucr.    try    rcr   measuri 
(EC  FT) 
FCr    L«l    TO    M 


;2*F3    i    r.3    IS    EW 


EISP    "I    air    wcrKinq     in    SCkT    rcr    ir.cc 


o       "    . 


;l 


2o  1 


F4(N)  =M(L)-G(L)/2+K*C  (L)/1G0 
F4$=vAL$(F4  (I  )  ) 
OUTPUT    717     ;E2$  ,F4$  ,ril$ 
wAIT    12*Tl+luG 

output  7o9   ;" vt3" 

a AIT     56 

ENTER    769     ;    A ( N ) 

NE  XT    i\ 

i-E£P 

GCLEAR     @    CLEAR 

SCALE    F 4  (  0 )  -  .  1*0  (L )  ,F4(1GG)  +  .1*0(L)  ,-.lf1.2 

XAXIS    0 ,G(L) /10  ,F4  (0)  ,  F4 (  100) 

YAXIS    F4 (0)  ,.1,6,1.2 

FGF    N=G    TC     100 

PENGP 

t-LC'l     F4  (N  )  ,A(N) 

NEXT    N 

MO  V  E    F  4  (  6 )  ,  -  .  1 

L A bEL    " MODE    " ; L ;  "     " ;  F  4  ( 0 )  ;  "    TC    "  ;  F 4  (  10 0) 

EiEP    "IF    YC0    '.-.ANT    A    CCFY    ENTER     1,    ELSE    ENTER    2"  ,@    E£ 

INPUT    Y 

IF    Y=2    THEN    9o7 

GCLEAR    @    CLEAR    @    SEEP 

EISP    "ENTER    FIGURE    #    TC    d£    PRINTED" 

INPUT    LS 

PLOTTER    IS     7  65 

P  EN    1 

SCALE    F4(0).-.2*O(L)  ,F4  ( 10  0)  +  .  1*0  (L)  ,-.3,1.3 

/.AXIS    U  ,C  (L)  /  iu  ,F4  (  0)  ,F4  ( 100) 

/AXIS    F4 (0)  ,  .1  ,o  ,1.  2 

PEN    2 

FCF    N=0    TC    106 

PEN  LP 
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962  FECI  F4(N)  ,A  (N) 
96  3  ME  XT  N 
9  64  t-Ei.LF  (o  PEN  1 
bob  EE1P  G, SIN (SO) 

960  FOE  X=*F4(0)  TO  F4(1CG)  STEP  C(L)/10  &  PENUP 
S  o  7  MO  v  E  X  ,  -  .  1 3 

968  LAEEL  V'AL$(X) 

96  9  NEXT  X 

970  LDIE  C  @  PENUP 

971  FOB  i=o  TO  1.2  STEP  .J  £  PENUP 
^72  i'iC  V  £  r  4  (  C  )  -  .  o  c  x  C  ( L  )  ,  Y 

97  3  LABEL  v'ALS(Y) 
9  7  4  NEXT  i 

975    LDIE     0     S    PENUP 

9  7  o    MG  v  E    F  4  (  G )  -h  .  1  *G  ( L )  ,-.21 

97/  LABEL  "FREQUENCY  IN  HS" 

978  LDIF-    0    Is?    FE}JuP 

979  MO v E    F  4 ( G ) - . 1*0  ( L )  ,  -  .  3 

scC  EAhEL  "PLOT  OF  RELATIVE  AMPLITUDE  vS  FREQUENCY,  MODE 

";L 

961  LDIF  G,SIN(9G.)  t  PENLP 
9b  2  MOVE  F4(  l).-.  17*0  (L)  ,  .  3 

S83  LAEEL  "RELATIVE  AMPLITUDE  " 

i>E4  LDIF  C  @  PENUP 

9c 5  MOVE  F4( U)+ .1*0 (L)  ,1.2 

bc6    LAEEL    "FIGURE    ";L£    @    PENUP    tr    PLCTTEF    IE    1 

ic7    CLEAR    lg    EEljt 

bco    DISP    "^c    ycu    want    c^   cnange    any  tnir.j, .    Ii    sc    enter    1, 

else    ^." 
^69    INPUT    it    CLEAR    e    I.EEP 
9SC    IF    \=2    THEN     120u 
IGoG    EISP    "Cnange    sensitivity?    Enter    1    icr    yes,    2    rcr    nc 

1G1C    INPUT    It    CLEAR    -    EcEP 

10  2  0    It     Y=2    THEN    1G6G 

1G3G    CLEAR    fc    EEEP 

lo-O   DISP    "Enter    new    sensitivity    setting" 

10  50    INPUT    A  5 

lOoO    CLEAR    fc    EEEP 

IG7G   EISP    "Cnange    rrecvvicth?    Enter    1    rcr   yes,    2    rcr    nc" 

108 G    IN P O.T    X (9    C EE A R    i§    EEEP 

1090    IF    Y  =  z:    THEN     1120 

11G0    EISP    "Enter    new    treqwiatn" 

1110    INPUT    O(L) 

1120    CLEAR    r-    i^EEP 

1130   DISP    "Cnange   center    rrec?    Enter    1    rcr    yes,    2    rcr    nc 
it 

1140    INPUT    Y&    CLEAR    (§    EEEP 

1150    IF    Y=2    THEN    llcG 

1 1 6  0   EISP    "Enter    new    center    tree" 
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no    c 


1170    IK'PLl    r  (L) 
llob    GCLEAF-,    ^    CLEAR 
119C   GCTC    771 
llyl    !    rrcn     120o    tc 
ter    irec   ar.u    tne  ^ 

12  U    CLEAF.    is    GCLEAE 

1201   D I S  F    "I    &n    w  cr  k  i  ng    c  n    i  r  e  c 

1210    H  ,bl  ,£4  fHl  ,H2  ,H3  ,H4*=0 

l^.^u    F  2  ,.u  ;  =  5u 

1 2 .:  0    F C I    l-  —  1    T C    1 0  0 

1240  IF    A(N)<A(£i)     TiiEK     12ttG 

1241  IF    A(tf)*A(h')     TEEK     1^2    £LFF     12 
12^2    rrlj'i     A  (r! ) 

12  50  A  6  { L )  *A  ( i\ ) 

126  0  F  b  ( I* }  *  F  4  ( 1\ ) 

12  7  0  tx~bi 

1 2  b  l  N  E  a  I    iNi 

1 2  S  0  A  7  ( L )  =A  o  (  L )  /  S  v  F  (  2 ) 


1680    is    tr,e    tirst    trv    a  l   qticir.c    cers 


r    ECRT   recce    ";L 


1310  IF  A7(L)=A(F)  TfiEK  1410 
L520  IF  A7(L)<A(Iv)  THEN  13  7l 

13  30  IF  A  (&)<£]  TEEN  l45o 


THEN  1 4 C U 


1340 

cl^(.) 

13  50 

a  i=r 

13  60 

GCTC    14  5  0 

13  7  0 

I  F    A  ( w  )  >  E  2 

138  0 

B2  =  A  (..) 

13  So 

ti  i.  =  i\ 

14  GO 

bLiC     14  50 

1410 

B  1  rE2«A  (  N  ) 

1^2  0  Kl rH2=N 


1440  GCTC  1,70 

14  50  uFa'i  N 

14  51  X=(A7(L)-£1)/(B2-£1) 

146 C  F7-X* (F4{H2) -E4   (HI ) )  +  F4 

14  7  0  FCF.  N=H  TC  100 

14 6 0  IF  A7(L)=A(F)  TflEK  1580 
L4 SO  IF  A 7  ( L )  >A  ( N )  THEN  15^0 
150  0  IF  A(M)>B3  THEN  1530 
1510  B3=A(N) 
1  5  2  u  H  3  =  i\ 
1530  GCTO  1620 

15  40  IF  A(N) <=3  4 
1550  £4  =  A  (N) 
156  0  H4=N 
15  7  0  GCTC  lo2C 
1560  F3  fE4«A(N) 
15b0  H 3 ,JB4-N 
lbOO  F  d  =  F  4  (  H  3  ) 


THEN  15  70 
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L610 

162o 
1621 

1630 
16*1 
I65u 
16b  o 
A  5 
16  7o 
16  71 
16  7  2 
1673 

16  cC 
169  C 
L6S  1 
169  S 
unc 
1700 
1701 
1702 
17G5 
17 10 

17  11 
1712 
1720 
1750 
17  55 
176  0 


GCTC  10  4  0 

NEXT  N 

X= ( A  7 (L) -£  4 ) / (B  3-6  4) 

Fb  =  -  (  X*  ( F  4  (H4  ).-r  A  (  H  .: )  )  )  +F  4  (fi  k) 

<*  (L)=F6  (L)/  (fa-t"7) 

PRINT  "MODE  ";L; 

PRINT  "CENTEF  ? 1 1^    IS  "  ;F6(L);"  AND  AMP  IS  ";P6lL)* 

PRINT  "C  IS  ";'*{L) 

CLEAR 

PRINT  USING  lb  7 3 

IMAGE  5/ 

NEXT  L 

CLEAF    6    GCLEAP 

CISF    "I    c.i    working    en    tiir.e    constant" 


!    Froir,    1700    t< 


7^0    the    largest   line    Constant    is   ic 


FCR    L=l   'IC   M 

T  (L)=^(L)/  (PI*f  6(L)  ) 

NEXT    L 

T=T(i) 

FCR    L=2    TO    M 

IF    T>T (L)    THEN     17  20 

T=T(L) 

ix  t^  1    L 

EEEF    &    CLEAR 

PR  IN  T    "  T  I  :•  E    CONST    >=    "  ;  T  *  10  C  0 

CISP    "Sec    ar.c    enter    new    tiire    constant    (in    fr.iilL-sec 


1" 

r 
i" 

i' 
l' 

i" 

i' 

l' 

i' 

l 

i' 

i' 

ianc( 

16CG 
ev.  , 
1601 
lbG2 
16  0  3 
16G4 
lcG5 
180o 


70 

cO 
9G 
91 
92 
S3 
&4 
95 

yo 

97 
96 
99 


DISP  "MUST  EE  GREATER  THAN  "  ;  T ;  "  sec" 

INPUT  '1 

Tl  =  * 

CLEAR 

PRINT  "NEW  TIME  CCNST  IS  ";T1 

PRINT  USING  1794 

image  4/ 

!  MEASURE 
I  MEASURE 
!  MEASURE 
!  MEASURE 

CISF  "I  err.  working  en  seance  re  cev  ,  4  cr  pes,  &  var 

a       " 

!    Frcrr.    lbOO    to    1924    is    tne    calc    cr    the    ncise,    st.    c 

anc    tne    norr.fcer    cr    cts.     usee. 

G,^=0 

F=(Fd  (2)  -co  (  1)  )/4+Fb{  1) 

F4$  =  vAL^  (F.) 

OUTPUT    717    ;H2$,F4$  ,H1$ 

WAIT    Tl*12+iC0 

OUTPUT    7o9     ;"vTl" 
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1607 
lcbc 
1809 
161u 
1611 
1612 
1813 
1814 
1815 
1816 
1317 
1618 

iais 

le  2b 

1821 

1  g  22 

1825 

18  3  b 

lo35 

1840 

io^5 

185u 

Lb55 

1080 

1861 

1864 

lo65 

icoo 

1667 

168  0 

lcol 

lcc5 

loo7 

18oo 

169C 

R£w" 

1691 

,5X, 

iti92 

"  (n.V 

1893 

1698 

1699 

19bb 

ure 

19C1 

19G2 

190  3 

190o 

19C7 


WAIT    1000 

FOP    ft  =  l    TC    10  0 

ENTER    7bS    ;    A(N) 

tv=W+A  (ft ) 

N EXT    ft 

»v=.v/  10  G 

FCF    ft  =  l    TC    10  0 

£=  (A  (N  )  -W  )  "2 

G=G+S 

NEXT    ft 

G*A5*Sv"R  (G/99  ) 

f  F I  ft  l    "  I  n  e   me  a  r.    is    "  ;  v. 

FRIft.l    "ine    stdr.carc    dev    is    "  ;  G 

PftIK'1    USING    lo22 

IftAGE    4/ 

FOR    L=]    TC    ft 

I ( L) =1 . 8  2*SwR (A6 (L) *A5/G) 

PRINT    "NO.    OF    POINTS    FGF    wCDE     ";L;"    IS    ";I(L) 

IF    I (L) <ibb    THE ft     lc5b 

I(L)=1G0 

IF    I ( L ) > 5 0    TH E ft    1  o  6  b 

I(L)=50 

I  (L) =IP(I  <L)  ) 

I (L}=24 

PRINT    "ftC.    CHOSEN    IS    ";!(£) 

PRINT 

ft  E  X 1    L 

CLEAR    is    E££P 

PFIft'i    USING    1681 

IMAGE     4/ 

PRINTER    IS    7b 1  ,1  32 

DIM    F  S.  L  2  C  ] 

Fi?="C-G    FAvIftE" 

PRINT  USING  16S1  ;  "  TIME"  ,  ,,;I£ftP  "  ,  "  BFES  "  ,  "ft  4","C.  F 
,  "AMB."  ,"Q"  ,"SKR"  f"C-AftP"  ,"*    PT"  ,F$ 

IftAGE  sA  ,5a  ,^A  ,6a  ,5A  ,  4X  ,  3A  ,6X  ,  "A  ,  SX  ,4A  ,  14X  ,  1A  ,9X  ,  3A 
5A  ,4X  ,4A  ,6X  ,10A 

PRINT  USING  lcS3  ;  "(SEC)","  J  C  )  "  , "  (HZ )  "  ,  "  ( Vrir.s)  "  , 
)" 

IftAGE     5A  ,4X  ,5A  ,  25X  ,  4A  ,  1GX  ,  fcA  ,3 OX  ,4A  ,/ 

PRINTER    IS    2 

CLEAR 

!  Frcir.  19CG  tc  496G  is  the  calculation  tor  the  meas 
inducing    center    fcrea    ana  Q 

15  =  1 

FOR    L=l    TG    ft 

DISP    "I    AM    IN    ftEASGFE    FOR    ftCCc     ";L 

F1=F6  (L)  -Fb  (L)  /C.  (L) 

F 2=F6 (L)+Fb (L)/w(L) 
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1906 

li>09 

1910 

1911 

IS  12 

1913 

1914 

1915 

1916 

1917 

1916 

1919 

1920 

1925 

1926 

19  27 

Ibla 

1929 

19  30 

1931 

19  32 

19  3  3 

19  34 

19  3  5 

19  3d 

19  37 

1936 

193  9 

1940 

1941 

1  94  2 
1943 
1950 
196b 
1970 
I960 
19  9o 
2iiuu 
2ulu 
2020 

2  o  30 
2o4U 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 


0=  (P2-P 1)/I (L) 

v  =  (  b  2 -iT  1)  /2 

ri=U 

IF  Ab(L) >. 3  THEN  1917 

Ad (L) =  2*^6 (L) 

Ao  (L) =2*At  iL) 

IF  AS  (L)  <A1  THEN  1911 

AS (L) =A1 

GCTC  19  2  5 

It  Ab(L) <.95  THEN  192  5 

AS (L) =Ac (L)/2 

Ao  (L)  =Ao  (L)  /2 

GCTC  19  17 

A$=VAL$(A6(L) ) 

Gui'PQT  717  ;A2$,A$,A3$ 

WAIT  200 

ENTEB  7  22  ;  16 

FCF  1^0    TC  I  (L) 

r4(:-i)=to(L)+u':i.-v 

F4$-VAL$ (F4 ( i ) ) 

H1$.="H2," 

H2$="rk" 

CLTPuI  717  ;H2$,F4$,Hl$ 

WAIT  12*11+10  0 

GbTPUI  709  ;"VT3" 

EN'IEF.  7  09  ;  A(N) 

NEAT  N 

EN  TEE  7  22  ;  15 

FC  K  N  =  1  'i  C  I  ( L ) 

IF  A (N) <A (H)  THEN 

IF  A(N) =A (ri)  THEN 

PRIN1  A(N) 

H=N 

N  E  X  T  N 

X=A  |H) 

X*(A(H  +  l).-A(H-l))/2 

2  =  (A  (H+l)  +A(H-1)  -2*X)/2 

F6(L) =~(Y/(2*Z)  ) 

A6  (L)=X  +  i*Fb  (L)+ii*Fb  (L)  ' 

Fb  (L )  =F4  (H)  4-F6  <L)  *  (F4  (Hi 

El  ,64  ,Hl  ,  H2  ,H3  rH  4=0 

82  ,B 3=5u 

A7  (L)  =Ab(L)  /SWE  (2) 

FGF  N«0  TC  H 

IF  A7(L)=A(N)  THEN  2^16 

IF  A7(L)<A(N)  THEN  2  3  70 

IF  A(N) <£]  THEN  2  45  0 

cl=A(N) 

B  1  =  N 

GCTC  24  5  0 


196u 

194  3  ELSE  1950 


-F  4  ( fi  - 1 J  ) 
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2  570 

'<.  3b0 

^3so 

24LO 

24  lu 

242u 
2430 

24  ^0 
2450 
2451 
246  I 
2470 

2<i6  0 

2^iG 

25  CO 
251C 
2520 

2  5  30 
2540 
255u 

25  6  0 
2570 
2560 
2590 

2olu 
2620 
2621 
2b  SO 
2b4u 
2b  7  3 
z.b  / 6 

26  77 
2678 
26  79 
2660 
2bo5 
26  6  9 
2690 
2  7  bU 
2710 
2711 
2712 
2720 
2730 
2740 
2  7  50 
2760 
2770 
2775 


IF  A(N) >E2 
£2=A  (iSl ) 
R2=N 


THEN 


CO 


GCTC 


.^50 


e  7  =  X  *  ( F  4  ( ti  2 )  -r  4  (hi)  )  +  £ 


EOF. 


2560 

25^u 


2620 

N  )  <  =  £4  THEN 


El ,E2=A(N) 

H  1  ,H2=*N 

F  7  =  F  4  ( H  2  ) 

GCTC  2470 

NEXT  N 

>.=  (A7  (L)  -Bl)/  (£2-131) 

(ill) 
I  (L) 
IF  A 7 (l) =A {.)     THEN 
IF  A 7  (£)  >A  (N)  THEN 
IF  A (N ) >£3  THEN  2  5 
B'3«A(N) 
H  3  =  N 
GCTC 
IF  A 
£4=A(N) 
E4=N 

GCTC  2o2  0 
£3  ,E4=A (N) 
hi  ,H4  =  N 
F5-F4 (n 3) 
GC'iC  2o4u 
N  E  X  'I  N 

X=(A7(L)-E4)/  (£3-54) 
FS=-(X*  (F4  (H4).-F4  (H  3)  )  )  +F4  (d  A. 
£  (L)  =ro  (L)/  (ta-t  7) 
PRINTER  IS  701  ,132 
CLEAR 


!  vAR¥  I* 

FCR  £  =  1  TC  2 

FCR  K=l  TC  3 

LISP  "I  AM  IN  RAVINE  FCR  w  OF  foCDE  ";L 

J  (1)=^(L) 

J  (2)=1.005*Q(L) 

J(3)«.S95*C{L) 

R=Ao (L) 

U=(1/0(K) )~2 

E (K) =0 

FCR  N»l  TC  I (L) 

C  (N) =R/  (J(K) *£vR(  (F4(W)/F6(L) -F b (L ) /F 4 ( N )  )  ~2  +  u)  ) 

D=(A(N)-C  (K)  )  "2 

E (K) =E  (K)+D 

NEXT  N 

CLEAR 
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2780 

261  3 
2614 
2676 
2677 
2876 
2679 
2660 
2661 
268  S 
2&SC 
2900 
2*10 
2911 
2*20 
2  9  30 
2940 


NEXT  K 
V=-( (E 

v(L)=u 
i 
i 
I 

I  VARY 
CLEAR 
FCR  K  = 
LISP  " 
j ( 1 ) =A 
J(2)=l 
J(3)=. 
U=(l/w 
£ (K) =U 
FOR  N= 
C  (N)=J 


(2)-E(3)  )/(2*(E(2)+E(3)-2*E(l)))) 

(  1)+V*.005*J  (1) 


ANP 


1    TO    3 

I    AM    IN    RAVINE    FCR   AMPLITUDE    OF    M5DE    " ; L 

6  (L) 

.0l2*A6 (L) 

998*Ab (L) 

(L))~2 

1    TO    I (L) 
(K)/(Q(L)*SCR((F4(N)/F6(L)-F6(L)/F4(N)) ~2+U) 


(A(N)-C  (N)  ) 

K)=E  (K)+C 

XT    N 

EAR 

XT    K 

-<<E(2).-E(3 

(L)=J  |1)+V* 


)  )/(2 

nn  -;* 


*  IF 


(E( 
(1) 


:)+£(3)-2*E(l)  )  )) 


2950  D= 
2960    E( 

29  70  NE 
2975  CL 
2  98  0  NE 
3020    V= 

30  30  A6 

3176  ! 

3177  I 

317  5  ! 

3179  !  VARY  FFEQ 
316  0  CLEAR 

318  1  FCR  K=l  TC  3 

31e9  DISF  "I  AM  IN  RAvINE  FOR  FREQUENCY  OF  I^.CLE  "  ;L 

3190  o (1)=F6 (L) 

3200  J (2)=.005*F6 (L)/Q(L)+F6(L) 

32  lu  J ( 3) =-( . U0  5*F6 (L)/Q(L) )+F6(L) 

3211  R=A6 (L) 

3212  0=  (I/O.  (L)  )  "2 
3220  E(K)=0 

32  30  FCR  N=l  TC  I(L) 

324  0  C(N)=R/(giL)  *SQR  (  (P4{N) /J(K).-J  (K)/F4  (N)  )  *2+U)  ) 

3250  D=(A (N) -C(N) ) "2 

3260  E (K) =E (K)+C 

3  27  0  NEXT  N 

3275  CLEAR 

3280  NEXT  K 

3  30  5  v=-(  (E(2).-E(3)  )  /  (  2*  (E  (  2)  +E  (  3)  - 

3 3 Co  F6(L)=0 (l)+V*.005*J(l)/g(L) 

3  307  .NEXT  E 

3306  ! 

3309  ! 


rE(l)  )  )  ) 
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3310 

3  311 

3  312 

3314 

3  3Io 

3  317 

3  3  IS 

3319 

3326 

3  321 

z   -  i.  <- 

z  :  *.  3 

3  3  26 

3  3  27 

3  32o 

335u 

3  3  o  0 

3  37  1 

3  3  co 

3  3  i  t 

34  6  u 

34o  1 

,  b'A  , 

3402 

I    LAST    SHOT    RAVINE     * 
CLEAR 

LISt    "I    AM    IN    LAST    SHGI    R" 
J  =A  b  ( L  ) 
u=(l/v(L) ) ~2 
£  =  U 

FOR    N=*l    TC    I  (L) 

C ( N ) =0 / (Q (L ) *SQR ( (F  4 (N) /F  6 ( L ) -F6 (L) /F4 ( N ) ) 
D=(A(N)-C  (N)  )  ~2 
E=S+D 
NEXI    N 
u  =J  *A5 
i=d/G 

u  I i'.   G$.  [  1 32  ] 
ASSIGN*     1    TC    "DATA" 
14=  TIME 
T7={T5+T6)/2 
t-  r  =  u 

PRINTff    1,15 
ASSIGN#     1    TO 


"2+U)  ) 


14  ,T7  ,P5,L  ,F6(L)  ,J  ,L(L)  ,Y  rA8(L)   , 

,4A    , 


G  •?  ~    ^L  ,  2  X  1 3D  .  .jD  i^X  (£  •  zL  l.  ,  ^.a  ,^.L  ,  •*  a  ,  .;LC  5  L  . 


3  L  .  4L  ,  3  a  ,  3  LC  3D  ,  3  >  ,4L  ,5a,  3L  ,  6  X  ,  D  .  4D  £  " 
PRINI    oSING    C$     ;    T4  ,1 1  ,B  5,L  ,F6  <L)   ,j,^(L)  ,i 


,AS  ( 


I5(E 

D  .  4DE 
L)  ,15 


3^1l    15=15+1 

3  ^  2  o    F  4  § = v  A  L  9  ( i?  6  ( L )  ) 

3  4  :  L    A$='vALv(At(L)) 

34  40    CuTPU'I     7i7     ;  A  2  5,  A  v,  A  3$ 

34 5 C    kvAIT    2 GO 

346  b  GOTPUT  7  17  ;K2$ ,F4$  ,fll$ 

3470  aAII  Tl* 16+luO 

B4BG  GOTPUT  7  6*  ;"vT3" 

34  90  ENTER  70S  ;  C8(L) 

4S79  PRINIER  IS  2 

CLEAR 

i%:h^'I  L 

PRINTER  IS  7b  1 

PRINT  OSING  4Sc4 

4  9  64  IMAGE  / 
4965  PRINTER  IS  2 
4*95  ! 

499b  I 
49  W  i 

4  99  6  I    TRACK 

5  u  U  6  T  9=  0 

5ulG  FCR  L=l  TC  M 

5011  DISP  "I  AM  TRACKING  MODE 


4*66 
49cL 
4962 
49  c  3 


5b 2  6  h 4£=VAL$ (r 6 (L) ) 
50  36  A$*VAL$ (Ac(L) ) 
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SO^b 

5C41 
5050 
5bbb 

5bol 
5u70 
50  6  0 
5090 
51bb 
511b 
5  i2u 
513b 
514b 
515o 
516b 
51/0 
521c 

^  i.  i.  Li 

52  3l 

524b 

52.5 

525o 

526  0 

5270 

529b 

5jut 

530  1 

5310 

4 

5  22u 


output  717   ;A2$,a$,ab$ 

aAI'I    20b 

CUTPLl    7  17    ; h2 v  ,1 ^ v  ,h 1^ 

SvAIl    1  1*16+ lb  b 

OUTPbU    7bb     ;"VT3" 

ENTER    7b 9    ;    A9 

IF    A9<=.9*Co(L)     THEN     5iib 

IF   A9>=l.l*Cb(L)    THEK     5110 

GCTC    52  4  5 

rib=l 

IF    A9> .  3    THEN    5  210 

Ad  (L) =2*A6 (L) 


AS =2* A 


bb 


o9 
899  9 

bObb 


5120 


I F    A  6  ( L )  <A  1    TBI 

A  8  (  L )  =A  1 

GCiC    5250 

IF    A9<.95    THEM    5  24  5 

A8 (L) =A6 (L) /2 

A9=A9/2 

GCTC    5  2  1b 

CL  EAR 

NEX1    L 

ENTER    7  22    ;    18 

IF    ALS (T8-T7) >=.2    THEN     19G2 

IF    T9<1    THEN     50  lb 

PRINTER    IS    2 

T4.=TIi<:E 

PPIN1  "GONE  TO  MEASURE  DUE  I'C 

GCTC     li>02 
CLEAR    @    GCLEAR 
DISP    "THE    END" 

£  i\  D 


A.iv.f  LT  IUDE 


TIME 


;  i 
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^PPENDrX  B 
A  sample  progran  is  lists!  for  obtaining  information 
from  the  data  taps  nd  utilizing  it  in  a  figurs  or  plot. 
The  key  spot  in  this  program  is  at  line  31,  the  reading  cf 
data  from  the  tape.  There  ace  a  maximum  of  S50  records  en 
the  tape.   In  this  case,  15  is  indicating  which  record  is  to 


be  accessed-     Each  of  these  rsccras   can  be  accessed   in 


random  manner.  Only  the  actual  value  is  stored,  not  the 
specific  variable  that  represents  it.  Thus  for  example,  if 
in  the  original  storage  program  the  values  F  =  1 ,  I  =  3, 
G  =  4,  and  3=7  were  stored,  Daly  the  actual  values  would 
be  stored  in  the  ori=r  given   (1,3r4,7).   tfhen  reading  these 


values  irom 


r  h  a  -  a  n  a 


ape,   the  1st  /ilue  might  be  called  a,   the 


second  Z,  etc.  The  important,  point  is  that  if  the  4th  entry 
in  the  record  is  dssired,  four  values  must  be  read  into 
memory  (read  n  values  to  get  2th  value)  and  whatever  the 
reading  program  calls  the  entry,  is  what  it  now  must  be 
referred  to  within  the  program. 


32 


1 

CLEAR  @  EEEP  @  GCLEAR 

2 

R' 

DISF  "HCw  MANY  POINTS  ARE 
■ 

3EI  KG  E  y/ALU 

AT 

EL?  ENTER 

INPU1  E 

4 

CLEAR 

^ 

CISF  "EN'lER  THE 

LOWEST  ' 

rEM 

P  DESIRED  C 

•  K 

FLCT" 

D 

INPUT  £1 

7 

CLEA^    @  EEEP 

b 

DISP  "EN'lER  THE 

HIGHEST 

TE 

:mp  c esi red 

ON 

PLC!" 

9 

INPUT  £2 

1U  CLEAR  £  BEEF 

1 

7  DISP  "ENTER  Th£ 

MODE  # 

DE 

:SIF£C" 

1 

6  INPUT  G 

1 

S  CLEAR 

2 

1  CISP  "ENTER  THE 

FIGURE 

# 

DESIRED  TC 

£E 

PRINTED" 

z. 

2  INPUT  £9 

4. 

2  CLEAR 

2 

7  DIM  Tl  (651)  ,.T2  ( 

851) 

2 

8  GCLEAR  <J  CLEAR 

NUMBE 


2  9   ASSIGN*     i    TC    "DATA" 

3  0    FOB    15=1    TC    E 

31    READ?     1,15     ;    1  4  ,T7  ,P  5  ,L  ,F  6 
22    IFL=G    THEN    4  0    ELSE    50 

40  Tl(  I5)  =  (F6/G)  ~2/10~6 

41  T2(I5) -T7+273.  16 

50  NEXT    15 

51  FLOTIER    IS    70  5 

52  PEN    1 

60  GCLEAR  @  CLEAR 

61  B 7 -FLOOR (£1+27  3.  lb) 

62  E6=CEIL (£2+273.16) 

70    SCALE    £7-10  ,£8+5  ,35  ,51 
60    XAXIS    40, 10, £7, £6 

50  YAXIS  £7 ,2 ,40 ,50 

51  LDIR  0, SIN (SO)  @  PENLP 

52  FCR  X=£7  TC  £8  STEP  10  @  PENUP 

53  MOVE  X  ,3  7. 5 

54  LAEEL  VAL$(X) 

55  NEXT  X 
96  LDIR  0 

57  FCR  Y=40  TC  50  STEP  2  @  PENUP 

58  MOVE  B7-4  ,Y 
9S  LABEL  VAL$(Y) 

100  NEXT  Y 

101  LDIR  0  (9  PENUP 

102  MOVE  B7+10 ,36 

1G3  LAEEL  "TEMPERATURE  (K)" 

104  LDIR  0,SIN(S0)  £  PENUP 

105  MOVE  £7-8  ,37 

106  LABEL  " (FREQ/M#) ~2  IN  KHZ~2" 
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107  LDIF    0    @    PEN  UP 

108  MOVE    £7-6  ,3  5 

1C9    LABEL    "PLCI    CF    fRE^~2    VS    TEMP ,MODE    "  ;  G 

110  LDIR    C    g    PEN  UP 

111  MOVE    B7  +  5 ,50 

112  L^EEL    "FIGURE    " ;ES 
14!?    PEN    2 

150  FOP   N=0    TO    B/ *-l 

151  P.«G+3*N 
15  2    PENUP 

152  PLOT    T2  (P)   ,T1  (P) 
154    NEXT    N 

15  5    PENUP 

IbO    CLEAri    @    GCLEAR 

lbl    REWIND 

io2    CISP    "riLL    CONE,     ITME    TO    GC    CN " 

170    EtE 
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